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Abstract 

The long-standing issues of determination of the mass distribution and nature 
of the center of our Galaxy could be probed by a lensing experiment capable of 
testing the spatial and velocity distributions of stars nearby and beyond it. We 
propose a lensing toy-model which could be a further evidence that a massive 
consensation {e.g. a neutrino condensation) is a good candidate to explain the 
data ruling out the presence of a supermassive black hole. 
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The puzzle to explain the nature of matter condensation at the center of our Galaxy 
is more than twenty years old problem |IJ. Various observational campaigns 0] have 
identified such a center with the supermassive compact dark object Sagittarius A* (Sgr 
A*) which is an extremely loud radio source. Detailed information comes from dynamics 
of stars moving in the gravitational field of such a central object. The statistical properties 
of spatial and kinematical distributions are of particular interest ||: Using them, it is 
possible to establish the mass and the size of the object which are (2.61 ± 0.76) x 10 6 M Q 
concentrated within a radius of 0.016 pc (about 30 Z<is) 0,[[5|]. 

More precisely, Ghez et al. [|J have made a campaign of observations where velocity 
measurements in the central arcsec 2 are extremely accurate. From this bulk of data, it 
is possible to state that a supermassive compact dark object is present at the center of 
Galaxy and, furthermore, it is revealed by the motion of stars moving within a projected 
distance of less than 0.01 pc from the radio source Sgr A* at projected velocities in excess 
of 1000 km/s. In other words, a high increase of velocity dispersion of the stars toward 
the dynamic center is revealed. Furthermore, a large and coherent counter-rotation, 
expecially of the early-type stars, is revealed, supporting their origin in a well-defined 
epoch of star formation. Besides, observations of stellar winds nearby Sgr A* give a mass 

accretion rate of — — = 6 x 10~ 6 M@yr 0. Hence, the dark mass must have a density 

~ 10 9 M & pc~ 3 or greater and a mass-to-luminosity ratio of at least 100M Q /L Q . The 
conclusion is that the central dark mass is statistically very significalt (~ 6 — 8a) and 
cannot be removed even if a highly anisotropic stellar velocity dispersion is assumed. 
Given that the majority of stars in a cluster are of solar mass, such a large density 
contrast excludes that the dark mass could be a cluster of almost 2 x 10 6 neutron stars or 
white dwarfs. As a first conclusion, several authors state that in the Galactic center there 
is either a single supermassive black hole or a very compact cluster of stellar size black 
holes H . The first hypothesis is supported by several authors since similar supermassive 
black holes have been inferred to explain the central dynamics of several galaxies as M87 
or NGC4258 |§. However, due to the above mentioned mass accretion rate, if 
Sgr A* is a supermassive black hole, its luminosity should be more than 10 erg s . On 
the contrary, observations give a bolometric luminosity of 10 37 erg s -1 . This discrepancy 
is the so-called "blackness problem" which has led to the notion of a "black hole on 
starvation" at the center of Galaxy. Besides, the most recent observations probe the 
gravitational potential at a radius larger than 4 x 10 4 Schwarzschild radii of a black hole 
of mass 2.6 x 1O 6 M |§ so that the supermassive black hole hypothesis at the center of 
Galaxy is far from being conclusive. 

On the other hand, stability criteria rule out the hypothesis of a very compact stellar 
cluster in Sgr A* ||. In fact, detailed calculations of evaporation and colision mechanisms 
give maximal lifetimes of the order of 10 8 years which are much shorter than the estimed 
age of the Galaxy |T0 |. 

Another viable and, in some sense more attractive alternative model for the super- 
massive compact object in the center of our Galaxy (and in the center of several other 
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galaxies) has been recently proposed by Viollier et al. [Hj. The main ingredient of the 
proposal is that the dark matter at the center of galaxies is made by nonbaryonic matter 
{e.g. massive neutrinos or gravitinos) which interacts gravitationally forming supermas- 
sive balls in which the degeneracy pressure of fermions balances their self-gravity. Such 
neutrino balls could have formed in the early epochs during a first-order gravitational 
phase transition and their dynamics could be reconciled with some adjustments to the 
Standard Model of Cosmology (for an exhaustive discussion of the problem, see |IT| ). 



Furthermore, several experiments are today running to search for neutrino oscillations. 
LSND |l!| finds evidence for oscillations in the v e — channel for pion decay at rest 



and in flight. On the contrary KARMEN [13] seems to be in contradiction with LSND 



evidence. CHORUS and NOMAD at CERN are just finished the phase of ('94-'95) data 
analysis. In any case, it is very likely that exact preditions for and — v T oscillations 
will be available at the end of millennium or in first years of the next. 

From all this bulk of data, and thanks to the fact that it is possible to give correct 
values for the masses to the quarks up, charm, and top, it is possible to infer reasonable 
values of mass for v e , v^, and v T . For our purposes, we are particularly interested in 
fermions which masses range between 10 and 25 keV/c 2 which cosmologically fall into 
the category of warm dark matterQ Choosing fermions like neutrinos or gravitinos in 
this mass range allows the formation of supermassive degenerate objects (from 1O 6 M 
to 1O 9 M ). As we said, the existence of such objects avoids to invoke the supermassive 
black hole hypothesis in the center of galaxies and quasars and it is able to justify the 
large amount of radio emission coming from such unseen objects. 

The theory of heavy neutrino condensates, bound by gravity, can be easily sketched 
T4| . Let us consider the Thomas-Fermi model for fermions. We can set the Fermi energy 



E F equal to the gravitational potential which binds the system, that is 

Ti 2 k 2 (r) 

' '"• - m v §{r) = E F = -m^(r ) , 



2m 



V 



where $(r) is the gravitational potential, hp is the Fermi wave number and $(r ) is a 
constant chosen to cancel the gravitational potential for vanishing neutrino density. The 
length ro is the estimed size of the halo. If we take into account a degenerate Fermi gas, 

1 /3 

we get kp{r) = (67r 2 n iy (r) / g u ) , where n„(r) is the neutrino number density and we 
are assuming that it is the same for neutrinos and antineutrinos within the halo. The 
number g v is the spin degeneracy factor. Immediately we see that the number density 
is a function of the gravitational potential, i.e. n v = /($), and the model is specified 
by it. If in the center of the neutrino condensate there is a baryonic star (which we 

1 A good estimated value for the mass of r-neutrino is 

m Vr = m v „ (—) =s UAkeV/c 2 , 



which well falls into the above range. 



2 



approximate as a point source), the gravitational potential will obey a Poisson equation 
where neutrinos (and antineutrinos) are the source term, i.e. 



A$ = -4nGm u n u . (2) 

Such an equation is valid everywhere except at the origin. We can assume, for the sake 
of simplicity, the spherical symmetry and define the variable u = r[$(r) — $(r )] then 
the Poisson equation reduces to the radial Lane-Emden differential equation 

fu _ ( AV2mjGg v \ 

dr 2 { 37rh 3 ) ^ ' 1 ' 

with polytropic index n = 3/2. This equation is equivalent to the Thomas-Fermi differen- 
tial equation of atomic physics, except for the minus sign that is due to the gravitational 
attraction of the neutrinos as opposed to the electrostatic repulsion between the elec- 
trons. If Mb is the mass of the baryonic star internal to the condensation, the natural 
boundary conditions are 

u(0) = GM B , u(r o ) = 0. (4) 
Recasting the problem in a dimensionless form, we have 

d 2 v v 3 ^ 2 
dx 2 -Jx ' 

and the boundary conditions 



(5) 



with the positions 



and 



Mr 

v (°) = j£> v(x o ) = 0, (6) 



u r 

v = rJwT x = n ' ( 7 ) 

GMq a 



3nh 3 \ 2/S /!7.2keV\ 8/3 



2.1377 g^Hjr. (8) 



4V2mig v G 3 / 2 MQ /2 J " \ rn v c 2 . 

We have to note that for m v = 17.2keV/c 2 the characteristic scale a, the corresponding 
of the Bohr radius for an electron bound to a nucleus, is, for a neutrino halo bound by 
a baryonic star, of the order of the average distance between stars. It strongly depends 
on neutrino mass. 

All the quantities characterizing the condensate can be written in terms of v and x 
(or u and r): 



n u ™>lg„ (2u\ 3 ' 2 n/x /6\ 2/3 vr 4 /^ 2 , x 



5/3 



(9) 



3 



which are, respectively, the gravitational potential, the number density, and the degener- 
acy pressure. As shown in [p] , the general solution of (H), or equivalently @, has scaling 



properties and it is able to reproduce the observations. In particular, it could well fit the 
observations toward the center of our Galaxy which estimate, considering the proper mo- 
tion (< 20 km sec -1 ) of the source Sgr A*, a massive object of M = (2.6 ± 0.7) x 10 6 M Q 
which dominates the gravitational potential in the inner (< 0.5pc) region of the bulge 
15 1 . In summary a degenerate neutrino star of mass M = 2.6 x 10 6 M Q , consisting of 



neutrino with mass m > 12.0 keV/c 2 for g v = 4, or m > 14.3 keV/c 2 for g v = 2, does 
not contradict the observations. Considering a standard accretion disk, the data are in 
agreement with the model if Sgr A* is a neutrino star with radius R = 30.3 Id (~ 10 5 
Schwarzschild radii) and mass M = 2.6 x 1O 6 M with a luminosity L ~ 10 37 erg sec -1 . 

Similar results hold also for the dark object (M ~ 3 x 1O 9 M ) inside the center of 
M87. 

Now the problem is: How much is the model consistent? Could it be improved at 
boundaries? Actually, due to the Thomas-Fermi theory, the model fails at the origin and, 
in any case, we have to consider the effect of the surrounding baryonic matter which, in 
some sense, have to give stability to the neutrino condensate. In fact, an exact solution 
of Eq.(§) or Eq.(||) is u(r) ~ r -3 from which $(r) ~ r~ 4 which is clearly unbounded from 
below. 

Let us now assume a thermodynamical phase where a constant neutrino number den- 
sity can be taken into consideration. This is quite natural for a Fermi gas at temperature 
T = 0. The Poisson equation is 

A$ = -AixGm u n u = cost . (10) 

which, for spherical symmetry, can be recast in a Lane-Emden form 

IT k 2 — Ucost, (11) 
r z ar \ ar J 

with polytropic index n = 0. The solution of such an equation is $(r) ~ r 2 , which 
is clearly bounded. For T / the above theory holds so that we get a solution of 
the form $(r) ~ r -4 . Matching the two results it is possible to confine the neutrino 
condensate. On the other hand, a similar result is recovered using the Newton Theorem 
for a spherically symmetric distribution of matter of radius R [JT6]] . In that case, the 
potential goes quadratically inside the sphere while it goes as $(r) ~ r -1 matching on 
the boundary. In our case, the situation is similar assuming the matching with a steeper 
potential. 

Assuming the existence of such a neutrino condensate in the center of Galaxy, it could 
act as a spherical lens for the stars behind so that their apparent velocities will be larger 
than in reality. Comparing this effects with the proper motion of the stars of the cluster 
near Sgr A* , exact determinations of the physical parameters of the neutrino ball could be 
possible. In this case, gravitational lensing, always used to investigate baryonic objects, 
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could result useful in order to detect a nonbaryonic compact object. Furthermore, since 
the astrophysical features of the object in Sgr A* are quite well known ||, accurate ob- 
servations by lensing could contribute to the exact determination of particle constituents 
which could be, for example, neutrinos or gravitinos. Besides, microlensing by cold dark 
matter particles and noncompact objects has been widely considered in literature ||17|| , 
being gravitational lensing independent of the nature and the physical state of deflecting 
mass. In fact, any gravitationally condensed structure can act, in principle, as a gravi- 
tational lens. Our heavy neutrino ball, being massive, extended and transparent, can be 
actually considered as a magnifying glass for stars moving behind it. If an observer is on 
Earth and he is looking at the center of our Galaxy (which is at a distance of 8.5 Kpc), 
he should appreciate a difference in the motion of stars since lensed stars and non-lensed 
stars should have different projected velocity distributions. In other words, depending 
on the line of sight (toward the ball or outside the ball) it should be possible to correct 
or not the projected velocities by a gravitational lensing contribution and try to explain 
the bimodal distribution actually observed 

Let us discuss the physical reasons why a heavy neutrino ball can be treated as a 
thick lens. 

For a static gravitational field, the refraction index is connected to the Newtonian 
potential by the equation 

n(r) = 1 - 2-^ (12) 



easily derived by 



2oo^l + 2-^; g ik ~-6 ik (l-2^\ ; (13) 



c 2 




assuming the weak field $/c << 1 and the slow motion approximation \v\ « c ||18|| . 
In this situation, almost all the usual geometric optics works. Our neutrino ball gives 
rise to a static gravitational field, it is an extended object and, due to neutrinos, it can 
be reasonably approximated by a "transparent" medium. In this case, it is essentially a 
thick spherical lens which can be replaced by two thin lenses at a distance "gT |19|| . It is 
easy to show, by elementary optics arguments, that this double dioptric system can be 
described by the equation 

11 1 d 

f ~ fi + h hh ' { } 

where are focal lengths. The relations with the gravitational field and the size of the 
neutrino ball are given by 

I = (n-l)i, £ = 1,2 (15) 
fi n 

and 7"j = R ~ d/2, where n (r) is the refraction index induced by the Newtonian potential 
R is the neutrino star radius. Since we are assuming a spherically symmetric distribution 
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of matter inside a radius R, the Newton theorem holds 



so that 



1 2 2 



with cij 2 = 4irGp/c 2 . The focal length is then given by [20] 



4Gp 



(16) 



(17) 



which has the same order of magnitude of / and fi$- If M — 2.6 x 1O 6 M and R ~ 30.3 Id, 
Lf oc ~ 70ly and we have to expect lensing effects approximately in this range for the 
stars behind the ball. Depending on the radius, Eq.([L2|) gives the refraction index so 
that the model is completely determined. Further considerations give also the range of 
validity of paraxial approximation. In fact, given a ray of light entering in the neutrino 
ball, it is possible to calculate, the angle of incidence a of the ray on the back surface, 
the angle of deflection S and the entering angle rj, (with respect to the normal) which 
produces the minimal deflection. See the Fig.l. 





Figure 1: sketch of the lensing effect by a neutrinos ball. 

By knowing these parameters, we could say if an intervening star undergoes this 
magnifying glass effect. Snell's law gives rasina = sin 7] so that the relation 

sin a = — sin rj < — , (18) 
n n 



holds. The angle of incidence a has to be smaller than the critical angle arcsin — so that 

n 

the incident light is reflected partially at the back spherical surface. The net effect is 
that we should lose a part of the luminosity of the star population behind the ball. As 

a = (rj — a) + x, we have 

5 = 7T - 2 X = 7T - 4a + 2ri; (19) 

which is the deflection angle. The minimal deflection is deduced in a straightforward 
way. We require 

^ = -4^ + 2 = 0, (20) 
drj drj 

, ... da 1,-,. ./1.x. i • , da cos ri 

which is — = -. Being a = arcsm(— sin 77), immediately we get — = and then 

dr\ 2 n dr\ ncosa 

1 1 • 2 4 2 
1 sin 77 = — cos 77 . 

n n 2 



Finally 



.2 " 2 



cos^ = — g— , (21) 

which relates the entering incidence angle with the refraction index. Summing up all 
these information, we should say if a given star undergoes or not a significant lensing 
effect behind the neutrino ball. 

Let us now take into account the projected positions (x, y) in the sky of the observed 
early- and late- type stars as reported in ||. The gravitational refraction index as given 
by Eq.(|l^) for a supermassive neutrino of mass M = 2.6 x 1O 6 M and radius R = 30.3 Id 
is n = 1 + 5 x 10~ 6 . From the above considerations, the angles rj and a are given by 

y sm 71 
rj = tan" 1 — , sina = (22) 

x n 

and then we can calculate the deflection angle 5 considering the refraction index given by 
the magnifying glass model (Eq.([H])) or given by the gravitational potential (p!2[). The 
mean values, taking into account the Genzel et al data || are 

So P t= 3.25 ± 0.89 arcsec , S gra v = 3.40 ± 1.91 arcsec, (23) 
for the sample of late-type stars, and 

5o P t= 3.48 ± 0.58 arcsec , S gr av= 3.49 ± 1.71 arcsec, (24) 

for the sample of early-type stars. In both cases the "optical" and "gravitational" results 
are in a good agreement. In other words, a magnifying glass model seems to reproduce 
the spatial distribution of stars behind Sgr A*. Fig. 2 shows such a spatial distribution 
for the late and early-type samples. The histograms for the deflection angle 5 are given 
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Spatial Distribution of Early- and Late-Type Stars 
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Figure 2: Spatial distribution of the Late-type and Early-type sample. 

in Fig.3. It is clear that a correlation exists between the deflection angle evaluated in 
both approaches. From gravitational lensing point of view, there is no relevant differ- 
ence between early and late type samples (see Fig.4). A more striking result concerns 
kinematics. The lensing effect of a possible neutrino ball at the center of our Galaxy will 
magnify stars up to 70lyrs behind, as we showed above, and their apparent velocities 
will be larger than in reality, as well (the effect is similar to that of looking at red fishes 
moving in a spherical water-jug of glass!). Taking into account the projected velocities, 
they will be corrected by lensing effects, i.e. 



v observed^ = y not lensed^ + ^ 



where 



Re-OeD*-^ c2 DoiDog , T - — t lensed (26) 

are respectively the Einstein radius and the time of crossing. Di s is the distance between 
lens and source; D Q i the distance between observer and lens; D os the distance between 
observer and source; Rb is the radius of the neutrinos ball and M its mass. We can 
assume 

D ls « D os « D os /2 « L foc = f . (27) 
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Early- Type Sample 

delta augje Willi Gaussian Fit 

Mean sd 

3,4373 0,88549 




Late -Type Sample 

delta angle v/iili c jnussian Fit 



Mean 



3.2213 0.5S11" 
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Figure 3: Deflection angle delta for Early- and Late type sample. 
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Fig. 5 shows y° bsert,ed (t) as measured by Genzel e£ a/. while the result of our evaluation 
is in Fig. 6. It is clear that if we take into account a lensing effect, the velocity dispersion 
{e.g. the sigma into the plots) of early-type stars becomes smaller and comparable with 
the late-type ones, considering that the standard deviation measured and reported in || 
is about 30 Km/ s. 

In conclusion, using such a magnifying neutrino ball, late and early-type stars could 
not have different spatial and kinematical distributions, the only difference should be if 
they are "behind" or "nearby" the neutrino ball from the observer point of view. In fact, 
as we have shown, early type-stars undergo a major lensing effect. If these considerations 
work, the central compact object could be investigated in such an alternative way and 
accurate kinematical and photometric data could give final answers on its size and nature. 



Acknowledgments 

The authors are grateful to G. Scarpetta and R.D. Viollier for fruitful suggestions and 
discussions which allow to improve the paper. 



9 



References 

[1] J. Oort, Ann. Rev. Astr. Ap. 15 (1977) 295. 

[2] R. Genzel and C.H. Townes, Ann. Rev. Astr. Ap. 25 (1987) 377. 

R. Genzel, D. Hollenbach and C.H. Townes, Rep. Prog. Phys. 57, 417 (1994). 

[3] K. Sellgreen, D.N.B. Hall, S.G. Kleinmann, and N.Z. Scoville, Ap. J. 317 (1987) 881. 

[4] A.M. Ghez, B.L. Klein, M. Morris, and E.E. Becklin, Ap. J. 509 (1998) 678. 

[5] R. Genzel, N. Thatte, A. Krabbe, H. Kroker and L.E. Tacconi-Garman, ApJ 472 (1996) 
153. 

[6] H.C. Ford et al. Ap. J. 435 (1994) L27. 
[7] Macchetto et al. Ap. J. 489 (1997) 579. 

[8] L.J. Greenhill, D.R. Jiang, J.M. Moran, M.J. Reid, K.Y. Lo, and M.J. Claussen, Ap. J. 
440 (1995) 619. 

[9] R.H. Sanders, Nature 359 (1992) 131. 

[10] E. Maoz, Ap. J. 447 (1995) L91. 
E. Maoz, Ap. J. 494 (1998) L181. 

[11] R.D. Viollier, Prog. Part. Nucl. Phys. 32 (1994) 51. 

N. Bilic, D. Tsiklauri, and R.D. Viollier, Prog. Part. Nucl. Phys. 40 (1998) 17. 

[12] D.H. White, Talk presented at Neutrino '98, Takayama, 4-9 June 1998. 

[13] B. Zeitnitz, Talk presented at Neutrino '98, Takayama, 4-9 June 1998. 

[14] R.D. Viollier, F.R. Leimgruber and D. Trautmann, Phys. Lett. B 297 (1992) 132. 
A. Merloni, R. Ruffini and V. Torroni, // Nuovo Cimento "5"113 (1998) 255. 

[15] A. Eckart and R. Genzel, Nature 383, 415 (1996). 

A. Eckart and R. Genzel, Mon. Not. R. Ast. Soc. 284, 576 (1997). 

[16] J. Binney, S. Tremaine Galactic Dynamics (Princeton Univ. Press.: Princeton 1987). 

[17] A.V. Gurevich, V. A. Sirota, and K.P. Zybin, Phys. Lett. A 207 (1995) 333. 
A.V. Gurevich and K.P. Zybin, Phys. Lett. A 208 (1995) 276. 
A.V. Gurevich, K.P. Zybin, and V.A. Sirota Phys. Lett. A 214 (1996) 232. 

[18] N. Straumann, General Relativity and Relativistic Astrophysics, Springer, Berlin, 1984. 

[19] P. Fleury and J. P. Mathieu, Images Optiques, Eyrolles, Paris, 1968. 

[20] S. Capozziello and G. Iovane, Grav. & Cosm. 5 (1999) 1. 



10 



Early -Type Sample 

Diff (Dellaopl - Della-grav) 




La le- Type Sample 

Diff (Dellaopl - Delia- a av) 



■1£ -1J3 43£ a.a 0.E i.n v 



delta vs. <lelta^yiavfLaie-Ty|>e Sample) 
delta_orav= 1.105 + 1.3683 1 delta 
Com elation: I = .89850 




3 

delta_opt 

delta vs. delta_yiavfEaily-Ty|>e Sample) 
delta jgrav = -1.559 + 1.4475 'delta 
Coiielation:i =.86711 




3 4 
delta_opt 



Figure 4: Correlation plots. 
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Figure 5: Histograms of observed velocities 
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Figure 6: Histograms evaluated velocities. 
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